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Abstract 

We report the results of a search for the 0*"*"^ pentaquark in the decay Q*~^^p where 

0*"'""'" — > pK~^ using 81 fb ~^ of data collected on the T(4S) resonance with the BABAR detector at 
PEP-II. We find an upper limit on the branching fraction of — > 0*"*"+^ where 0*+"*" — > pK^ to 
be 1.5 X lO"'^ for 1.43 < m{@*++) < l.SSGeV/c^, 2.4 x lO"'^ for 1.85 < m(e*++) < 2.00GeV/c^ 
and 3.3 x lO"'^ for 2.00 < m(e*++) < 2.36GeV/c^ at 90% confidence level. Ah results are 
preliminary. 
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1 INTRODUCTION 



Recently several experimental groups have reported observations of a new, manifestly exotic baryon 
resonance, called the 0'''(154O) [P, with an unusually narrow width (r < 8MeV/c ). These results 
have prompted a surge of pentaquark searches in experimental data of many kinds j2j. In this paper 
we will concentrate on the exclusive search for pentaquarks in the decay of B mesons. Following 
the observation of the decay B'^ ppK~^^ [SI 15 it was suggested that this decay might include 
events of the form B~^ — > Q*~^^p where 0*"'"+ is an I = 1, I3 = 1 pentaquark [3]. O*"'""'" would 
be a member of the baryon 27-plet with quark content uuuds. It has been predicted to lie in the 
region 1.43 — 1.70GeV/c^ in the pK^ invariant mass of B^ ppK^ candidates and to have a 
width of 37 — 80 MeV jHl • The pK^ cross section is nearly purely elastic in the region of interest 
so a resonance would follow the Breit-Wigner form, with a peak cross section of about 25 mb if the 
resonance is at 1.7GeV/c^ and even larger if the mass is less. The cross section is measured to be 
about 12 mb at center-of-mass energies spaced by about 15 MeV so its width would need to be 
considerably less than 15 MeV to have escaped detection. Our limits will not depend on the 0*"'""'" 
width in any significant fashion. We will search for 0*"*"^ in the mass region up to 2.36GeV/c^. 

2 THE BABAR DETECTOR AND DATASET 

We use data collected on the T(4S) resonance with the BABAR detector at PEP-II to search for 
Q*++_ 'pj^g data sample contains 89 million BB pairs, corresponding to an integrated luminosity 
of 81fb~^ on the T(4S) resonance. An additional 9fb~^ of data, collected 40 MeV below the 
resonance peak (referred to as off-peak data), are used to study the background from light-quark 
and cc production. 

A detailed description of the BABAR detector can be found elsewhere [H]; only detector compo- 
nents relevant to this analysis are mentioned here. Charged-particle trajectories are measured by 
a five-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift-chamber (DCH), oper- 
ating in the magnetic field of a 1.5-T solenoid. Charged particles are identified by combining the 
measurements of ionization energy loss {dE/dx) in the DCH and SVT with angular information 
from a detector of internally reflected Cherenkov light (DIRC). Photons are identified as isolated 
electromagnetic showers in a CsI(Tl) electromagnetic calorimeter. 

3 ANALYSIS METHOD 

We require that charged tracks have a minimum transverse momentum {pt) of 0.1 GeV/c, at least 
12 hits in the DCH, and that they originate from the interaction region point within 10 cm along 
the beam direction and 1.5 cm in the transverse plane. 

The kaon and proton particle identification is based on dE/dx information from the DCH and 
SVT for Pt < 0.7GeV/c or the measured Cherenkov angle and the number of photons observed in 
the DIRC for pr > 0.7GeV/c. 

The B candidate is formed from the proton, the anti-proton and the kaon candidates. Two 
kinematic variables are used to isolate the B~^ — > ppK~^ signal taking advantage of the kinematic 
constraints of B mesons produced at the T{AS). The first is the beam-energy-substituted mass, 
mES = [(-^cA/j/^ + Pi ''PbY l^i ~'P\Y^'^ ■• where Ecm is the total center-of-mass energy of the e+e" 
^Charge-conjugates are assumed throughout the paper. 
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collision. Here, the four-momentum of the initial e^e~ system is {Ei, pj and is the momentum 
of the reconstructed B candidate, both measured in the laboratory frame. The second variable is 
AE = E*g — Ecm/'^^ where E'^ is the i?-candidate energy in the center-of-mass frame. 

Several topological variables provide discrimination between the large continuum background 
{e^e~ qq, where q = u,d, s, c), which tends to be collimated along the original quark direction, 
and more spherical BB events. In order to suppress the dominant continuum background we use a 
linear combination (a Fisher discriminant) of the following four event-shape variables: cosO^^., the 
cosine of the angle between the thrust axis of the reconstructed B and the beam axis in the center- 
of-mass frame; cos^^^^, the cosine of the angle between the momentum of the reconstructed B 
and the beam axis in the center-of-mass frame; and zeroth- and second-order Legendre polynomial 
momentum moments, Lq = J2i \Pi \ and L2 = J2i |Pi |[(3cos^ Othrs i ~ where p* are the center- 

of-mass momenta for the tracks and neutral clusters that are not associated with the B candidate, 
and Othrg i are the angles between p* and the thrust axis of the B candidate. We optimize the 
Fisher discriminant coefficients for the best background and B^ ppK^ signal separation using 
off-resonance data and B^ ppK^ simulated events that are distributed uniformly in phase- 
space {B~^ — > ppK'^ signal Monte Carlo). These event topology requirements retain 67% of B~^ — s- 
ppK^ signal while removing 94% of continuum background. We expect 94% of the combinatoric 
background to come from continuum events and the remaining 6% from BB events. 

The Fisher discriminant, niEs, and AE cuts are optimized to maximize the statistical sensitivity 
of the B'^ ppK~^ signal, defined as S/^J S + B, with S and B being estimated numbers of 
i?+ — > ppK^ signal and background yields in the Monte Carlo simulation respectively. We assumed 
the B^ ppK^ signal branching fraction of (5.66!to'57 ± 0.62) x 10^^ in the optimization. The 
B+ ppK+ signal region is defined to be 5.276 < ruES < 5.286 GeV/c^ and \AE\ < 0.029 GeV 
(signal box) and the tues sideband region is taken to be 5.20 < ruES < 5.26GeV/c^ and \ AE\ < 
0.029 GeV for the combinatoric background studies. 

The main source of BB backgrounds is the b — > ccs transitions, where B'^ — > XccK'^ , Xcc — > pp 
and Xcc = Vc, J/'fp, tpC^S), Xc0,i,2 (so-called "charmonium background"). We expect 72 it 10 events 
of this type in the signal box region. To check for additional BB backgrounds that might peak in 
the B^ — > ppK^ signal region, we study generic BB Monte Carlo as well as a set of samples of 
exclusive B decay simulated events for potential charmoniumless backgrounds. The expected BB 
"charmoniumless" background contribution is less than one event in the signal box region. 

4 SYSTEMATIC STUDIES 

Systematic uncertainties in the analysis are described below and are summarized in Tabled 

The B'^ @*^'^{pK'^)p signal efficiency is computed with B^ — > ppK^ simulated events, 
reconstructed and selected using the same procedure as for the data. We apply small corrections 
determined from data to the efficiency calculation to account for the overestimation of the tracking 
and particle-identification systems performance. The resulting efficiency, as a function of mpx+, 
is shown in Fig. ^ A systematic uncertainty is assigned to each correction to account for the 
limited size and purity of the control sample used in computing that correction. For example, for 
the kaon identification, we correct the simulation using a pure sample of D*~^ — > 7r~^D^ decays with 
— > K~TT~^ and for the proton identification we use a sample of yl — > pTT~ . Conservatively, we 
take the size of the correction applied as our systematic error. 

In addition, after all the corrections, we compare our B^ ppK~^ signal simulation to a 
control sample with similar kinematics and final state topology {B^ — > J /'il){e^ e~)K~^)^ in order 
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Table 1: Systematic uncertainties for the branching fraction of Q*^^{pK~^)p without(with) 

background subtraction. 



Type 


% BF 


i?-counting 


1.1 


Tracking 


2.4 


PID 


6.0 


Event Shape 


2.0 


Signal Box Cut 


2.5 


Monte Carlo Statistics 


1.1 


Background subtraction 


(1.1) 


Total 


7.4(7.5) 



to quantify the ability of the simulation to model the kinematic and event-shape variables used in 
the event selection. The small residual differences in the efficiencies at the cut value are assigned 
as systematic uncertainties affecting the selection procedure. 

For the calculation of the systematic error due to the background subtraction we decrease the 
i?-backgrounds by the uncertainties in their branching fractions [7j. The change in the upper limit 
for the new background estimation is 1.1% which we take as a systematic error. 

The systematic error also comprises the uncertainties from the determination of the number of 
BB pairs. We assume that the branching fraction of the T(45) into BB is 100%, with an equal 
admixture of charged and neutral B final states. We do not include any additional uncertainty due 
to these assumptions. 

5 RESULTS 

The distribution of events in the ruES — plane is shown in Fig. [3 We see 212 events in 
the signal box region. To extract the number of B'^ ppK~^ signal events we loosen signal box 
cuts on AE and fit the AE projection for 5.276 < ruES < 5.286 GeV/c^ and |A^| < 400 MeV 
with a single-Gaussian distribution for B^ — > ppK~^ signal and a first-order polynomial for the 
background. From that fit we estimate that the 212 total events comprise 40ib2 combinatoric 
background events and 188ibl7 B^ ppK'^ signal events, including 68ibl0 events originating 
from charmonium decays to pp. 

The Dalitz plots for the events in the signal box (212 events) and the sideband region (368 
events) are shown in Fig. |51 Note that the relative phase-space, or the fraction of Fig. |31^bottom) 
in the signal box, is 0.104. The distributions in Fig. 01 are not efficiency-corrected. The Dalitz plot 
for the events in the signal box (Fig. ^ shows a threshold enhancement in the pp mass spectrum, 
as well as three clear bands corresponding to r/c, J/V' and i^i'^S) events. The background events 
tend to lie on the edges of the Dalitz plot because they are dominated by inclusion of random soft 
tracks. 

As we are interested only in the low mpj^+ region the following figures will be limited to rnpx+ 
up to 3.4GeV/c^ or the total of 75 events in the signal box region. It is convenient to represent 
data in two different ways: in Fig. |^ we separate the events into those inside the charmonium 
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Figure 1: The — > ppK^ signal reconstruction efl&ciency(left) and the detector resolution(right) 
as functions of mpx+ 
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Figure 2: ruES ~ distribution of on-peak data reconstructed in the ppK~^ mode. The 

small box (blue) is "signal" box: 5.276 < tues < 5.286 GeV/c^ and \ AE\ < 29MeV; and the large 
box (red) is "sideband": 5.20 < niEs < 5.26GeV/c2 and \AE\ < 29MeV. 
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Figure 3: Dalitz plot of on-peak data reconstructed in ppK^ mode. Events in the signal box 

region (top), sideband region (bottom). Note that these distributions are not efficiency-corrected. 
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Figure 4: The mpx+ distribution for data events in ppK~^ signal box: events in the char- 

monium region 2.85 < mpp < 3.15GeV/c^ (sohd), events outside the charmonium region (dashed). 
Note that these distributions are not efficiency-corrected. 
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Figure 5: The mpx+ distributions for data reconstructed as B'^ — > ppK^ (dots), niEs sideband 
(empty squares) and, Monte Carlo, B^ — > 7]c{pp) {dashed line) and B^ J /'>p{pp)K~^ (solid 
line). Note that these distributions are not efficiency-corrected. 
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Figure 6: Upper Limit on the branching fraction of Q*^^{pK^)p at 90% confidence level 

with the assumption of no background(dashed), with background as determined from data and 
Monte Carlo (solid). The systematic error correction is included in the limits. 

window and those outside, where as in Fig. [SJwe emphasize the different background contributions 
to the data. 

We search for B*^^ pentaquark in the pK^ mass spectrum, shown in Fig. [l] The binning 
corresponds to 4- crpx+, where o"px+ is the detector resolution shown in Fig. fright). The average 
B+ Q*++{pK+)p signal efficiency is (17.0 ±0.2)% for 1.43 < m,pK+ < 2.40GeV/c2. We observe 
no events for mpx+ < 1.85GeV/c^. 

The background contributions are shown in Fig. |S1 The rnpx+ distribution of the combinatoric 
background is obtained from the events in the data niEs sideband region and is scaled to the 
expected number of the combinatoric background events in the signal box. The shape and amount 
of — > rjc{pp)K~^ and J / tp{pp)K'^ background contributions are determined from the 

simulation and scaled by their respective branching fractions 7 . There is no contribution to the 
background from r]c{pp)K~^ for mpx+ < 1.80GeV/c^ and B^ J/ip{pp)K~^ for mpx+ < 

1.75 GeV/c2. 

To set an upper limit at 90% confidence level on the branching fraction of B~^ — > Q*^^ {pK~^)p 
we count events in each of the rnp]^+ mass bins in Fig. |1] assuming that all the events observed are 
B~^ — > Q*^~^ {pK~^)p signal events. To simplify the presentation of the upper limit on the branching 
fraction as a function mpx+ we assume that number of events in each of the bins in mpx+ is equal 
to the maximum number of events per bin for each of the mpj^+ regions (see Table |2) ■ 

We use two methods to determine the upper limit. In the first one we assume that there 
is no background contribution. We calculate from Table 31.3 [7j the Bayesian upper limit for 
90% confidence level as a function of mpx+ assuming Poisson-distributed events in the absence 
of background. The resulting values are shown in Fig. El and given in Table [21 To account for 
systematic errors we increase the upper limit by the total systematic error (7.5%). 
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Table 2: The upper limit for the branching fraction of Q*^^{pK^)p as a function mpj^+ 

without(with) background subtraction. 



Mass Region, 


Maximum Events Observed 


BF UL (lO-'O @ 


BF UL (lO-'O @ 


GeV/c2 


in any mpx+ bin 


90% CL without bkg 


90% CL with bkg 


1.43< mpj^+ <1.85 





1.63 


1.49 


1.85< mpK+ <2.00 


1 


2.76 


2.40 


2.00< mpK+ <2.36 


2 


3.78 


3.28 



To calculate the upper limit in the presence of background we use a tool described in 0. It 
uses toy Monte Carlo technique to calculate an upper limit in presence of uncertainties on the 
efficiency and the number of expected background events. We assume all the systematic errors but 
the systematics on background and i?-counting to contribute to the uncertainty on the efficiency 
(7.3%). To estimate the number of expected background events we fit a first-order polynomial 
to the pK~^ mass spectrum of the combinatoric background events as well as r]c{pp)K^ 
and — > J/'ip{pp)K^ Monte Carlo events (so-called peaking 5-background) . The uncertainty 
on the background comes from the statistical error on the fit as well as the systematic error on 
the background. The resulting values of the upper limit as a function of mpj^+ increased by the 
systematic error on i?-counting (1.1%) are given in Table I^J and shown in Fig. El 

6 SUMMARY 

Using 81 fb^^ of on-peak data accumulated by the BABAR detector, we set an upper limit at 
90% confidence level on the branching fraction of B'^ Q*~^^{pK^)p to be 1.49 x 10~^ for 
1.43 < m(G*++) < 1.85GeV/c2, 2.40 x 10"^ for 1.85 < m(e*++) < 2.00GeV/c^ and 3.28 x lO"'^ 
for 2.00 < m(e*++) < 2.36GeV/cl 
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